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We introduce an infectious default and recovery model for N obligors. The obligors are 
assumed to be exchangeable and their states are described by N Bernoulli-type random 
variables S%{i = 1, ■ • • , N). They are expressed by multiplying independent Bernoulli vari- 
ables Xi,Yij and Yy, and the default and recovery infections are described by and Y^. 
We obtain the default probability function P(k) for k defaults. By considering a contin- 
uous limit, we find two nontrivial probability distributions with a reflection symmetry of 
Si *-> 1 — Si. Their profiles are singular and oscillating and we theoretically investigate it. 
We also compare P{k) with an implied default distribution function inferred from the quotes 
of iTraxx-CJ, which is a portfolio credit derivative of Japanese 50 companies. In order to 
explain the behavior of the implied distribution, the recovery effect may be necessary. 
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1. Introduction 

The cooperative phenomena, especially phase transitions, have been extensively studied 
and continue to be important subjects until today. They have provided universal paradigm 
for physics, sociology, and economy. The economical systems composed of a large number of 
interacting units have been studied from this viewpoint. 1 ' 2 Recently, systemic failure problems 
are being focused upon in econophysics, 3-7 financial engineering, 8-12 and computer engineer- 
ing, 13 and many probabilistic models have been proposed. This has been motivated by the 
fact that the description of systemic failures is necessary to control and manage them. An- 
other motivation is that credit risk markets are now growing, and the pricing of products is 
an immediate concern. 14 For this purpose, it is necessary to develop probabilistic models that 
can describe credit risks. 

The difficulty in the description of systemic failures arises from the fact that they are 
not independent events. If they are independent, the description is very easy and we only 
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need Bernoulli-type random variables Si denoting the element i's failure or not by Si = 1 or 
Si = respectively. However, there are many phenomena wherein the "correlation" between 
the failure events is very important. For example, in a network of storage systems, if a node 
fails, the failure can propagate to other nodes. 

In credit risk markets, the same type of risk propagation is found to occur. A percolation- 
type probabilistic model was proposed to describe bank bankruptcies, where interbank de- 
posits lead to collective credit risks. The probability of I failures obeyed the power law 
P(l) ~ / 1_T with the Fisher exponent r at its critical point. In a study, 8 a default infec- 
tion mechanism was proposed to describe the risk-dependency structures. The constituents 
are obligors, and the risk is whether he (or she) can refund before the expiry date. Such a risk 
is called a default risk. Davis and Lo introduced independent Bernoulli- type random variables 
Yij, which describes the infection from a bad obligor j to a good one i. They explicitly ob- 
tained the probability function for k defaults, P(k). They estimated the effect of the default 
correlation on P(k). 

One of the crucial problems with these studies is that their it is difficult to describe 
whether P(k) values do describe the empirical default distribution P(k) or not. Because of 
the relative scarcity of good data on credit events, it was impossible to compare the models 
by using the empirical data. Recently, from the market quotes on credit risk products, it 
becomes possible to infer the default distribution function. 12,15 We can compare and calibrate 
the probabilistic models. In the present paper, we generalize the model proposed by Davis 
and Lo by introducing a recovery effect. We compare the default distribution P(k) with an 
implied value of that of the credit market and calibrate the model parameters. With regard 
to the bulk shape, we see that the calibrated P(k) value looks similar. 

The outline of the present paper is as follows. Section 2 gives a brief introduction to 
the infectious default model proposed by Davis and Lo, 8 and we modify it by introducing a 
recovery process. We obtain the default probability function P(k) for k defaults. In section 3, 
we take the continuous limit of P(k) with finite Pd and non zero correlation p > 0. We find 
two non trivial probability distribution functions with a reflection symmetry. They exhibit 
oscillating behaviors and we investigate the mechanism. We compare the model distribution 
function P(k) with the market implied function in section 4. Section 5 summarizes our results 
and future problems are discussed. 

2. Infectious default model 

We consider N exchangeable obligors whose states are described by random variables 
Si(i = 1, 2, • • • , N) such that Si = 1 if obligor i defaults and Si = otherwise. Here, the term 
"exchangeable" means the non-dependency of the joint probabilities P(S\, S2, • • • , Sn) on the 
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exchange of Si <-> Sj for any pair of The number of defaults is 

K = S 1 + S 2 + --- + S N . (2.1) 

The value of Si is determined as follows. For % = 1, • • • ,N and j = 1, • • • , JV with j ^ i, let 
Xj, y,- be independent Bernoulli- type random variables with probability function 

Prob.pQ = l] = p, 

Prob.[Fy = l] = q. (2.2) 

5i are defined as 

s i = x i + (i- Xi)(i - n^(i - y,x,)). (2.3) 

Here, Xi is the internal state variable which describes whether the obligor is in a good state 
(Xi = 0) or not (Xi = 1). Si is also the state variable that describes whether the obligor is 
defaulted (Si = 1) or not (Si = 0), this is determined by not only the internal state but also 
the external environment. If Xi = 1, Si = 1 and the obligor is defaulted. Even if Xi = 0, 
obligor i can be defaulted. represents the influence of another bad obligor (Xj = 1) on 
obligor i. A default infection from a bad obligor j takes place if y,- = 1 and Xj = 1. Sj 
becomes 1 and the obligor is defaulted. The second term of eq.(2.3) represents this effect. 

We introduce a supporting effect from other good obligors in addition to the default 
infection. In fact, it may occur that a good obligor supports other bad obligors and the 
latter can circumvent their defaults. We introduce new independent Bernoulli-type random 
variables Y-j in addition to eq.(2.2). For i = 1, ■ ■ ■ ,N and j = 1, • • • ,N with j / i, they have 
the probability function 

Prob.[y/. = l]=</. (2.4) 
We introduce the following model equation for Si'. 

Si = XiU^i - y/.(i - Xj)) + (i - Xi)(i - n jyi (i - y,x,)). (2.5) 

Eq.(2.5) reveals that even when Xi = 1, if Xj = and Y(- = 1, obligor i is supported by obligor 
j and avoids being defaulted. We note that eq.(2.5) has a default, non-default symmetry. We 
get 1 — Si by substituting X j — ► 1 — Xi and Yij <-> y/.. This model can be reduced to the 
original infectious default model by substituting Y-j = into eq.(2.5). 
The probability distribution function P(k) for k defaults is given by 

k N-k 

P(k) = Prob.[K = k] = N C k x^E °Wk( l >rn) (2.6) 

1=0 m=0 

where 

a 7V,fc l'' m J — k^l x N-k^m X P {L-P) 

x (l — q') l ( k + m - 1 ) (I — g^m(N-k-m+l) 

x (l - (1 - qjN-k-m+ljk-l^ _ (1 _ q^k+m-l^N-ik+m) _ ^ 
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Fig. 1. Pictorial representation from internal states {Xj} to a conclusive state of k defaults and N — k 
non-defaults. 



We explain the derivation of eq.(2.7). In Fig.l, there are ./V obligors, k obligors are defaulted 
and N—k obligors are non-defaulted. The k defaulted obligors are classified into two categories: 

(A) and (B). (A) contains I bad obligors, which are never supported by other good obligors. 

(B) contains k — I good obligors, which are infected by other bad obligors and thus get 
defaulted. The number of different possible combinations of / items from k different items is 
kCi. Further, there are two categories (C) and (D) for the N — k non-defaulted obligors. (C) 
contains N — k — m bad obligors and (D) contains m good obligors. The N — k — m bad obligors 
are supported by other good obligors and they are prevented from being defaulted. The m 
good obligors are never infected to be defaulted. The number of different possible combinations 
of m items from N — k different items is jv-fcCm- In other words, the conclusive k defaults 
and N — k non-defaults are made from the N — k — m + I bad obligors and k — I + m good 
obligors in the internal configuration by the infection and recovery mechanism. The internal 
configuration is realized with probability p N - k - m + l (\ _ p^k—l+m_ i DaQ i obligors among the 
N — k — m + 1 obligors are not supported by the k — I + m good obligors; this probability is 
given by {\ — q'^ k ~ l+m ) 1 . m good obligors are never infected by N — k — m + l bad obligors; this 
probability is given by (i_g) m ( Ar - fc - m +0_ k, — \ good obligors must be infected by N — k — m+l 
bad obligors; this probability is given by (1 — (1 — qr) m-i-Z^A: z _ ^ _ m DaQ i obligors must 
be supported by k — l + m good obligors; this probability is given by (1 — (1 — q'^ k + m - l ^ N - k - m _ 
Therefore, the probability of k defaults and N — k non-defaults from a configuration (l,m) is 
given by oP^ q ^ (l,m), as shown in eq.(2.7). We obtain P(k) as the summation of (I, m) 
over I, m. 
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The expected value of the number of defaults K is 

< K >= N[p(l - q'il-p))*- 1 + (1 - p)(l - (1 - qpf- 1 )}, (2.8) 
and the default probability Pd is given as 

P d =< K > /N = p(l - q>(l - p)) N ~ l + (1 - p)(l - (1 - qpf- 1 ). (2.9) 
The variance is 

a\ =< K > +N(N - l)f3 p r f' q '- < K > 2 , (2.10) 



where 



= p 2 {l-2q'(l-p)+q' 2 (l-pf} N - 2 

+ 2p(l-p)(l-q'){(l-q'(l-p)) N - 2 
-(l- q )(l- q '(l-p)-pq) N ^} 
+ (1 - pf[l - 2(1 - pq) N ~ 2 + (1 - 2 M + pg 2 )^- 2 ] (2.11) 
and the correlation coefficient is given by 

Pn_ £k. (2.12) 

^ Pd(l-Pd) K ' 

We find that there are multiple solutions (p,q,q') corresponding to a value of P^. In 
particular, for large N, there are three solutions. For example, there are three solutions p = 
0.808310,0.5, and 0.191680 for N = W0,P d = 0.5, and q = q' = 0.05. On the other hand, 
there is only one solution p = 0.5 for N = 50, P d = 0.5, and q = q' = 0.05. This is because 
for arbitrary q, q' ^ 0, Pd behaves as that in Fig. 2. Pd(p, q, q') starts from at p = to 1 at 
p = 1. For intermediate values of p, Pd rapidly increases to 1 and then decreases to near 
p = 1 in the large N limit. Thereafter, Pd again increase rapidly to 1 toward p = 1. Such 
a behavior can be explained by eq.(2.9). There are three p solutions corresponding to a Pd 
value; they are referred to as left, middle and right solutions according to the order of p. The 
parameter region (q, q') in which there are three solutions of p expands with N; further, for 
the limit TV — ► oo, it covers the entire parameter space (q > 0, q' > 0). 

The profiles of the three solutions are shown in Fig. 3. We set N = 50, Pd = 0.5, and 
q = q' = 0.2. The three solutions are realized at p = 0.079281 (left), p = 0.5 (middle), and 
p = 0.920719 (right). The profiles of the probability distribution functions of the left and 
right solutions are reflection symmetric. The origin of the symmetry arises from the reflection 
symmetry of eq.(2.5). P(k) for the middle solution (p = 0.5) has a symmetrical profile and is 
almost a binomial distribution Bi(50, 0.5). 
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Fig. 2. Plot of P d vs. p. We set N = 50, 100, Fig. 3. Plot of P(k) for P d = 0.5 and q = 
and 500, q = 0.3, and g' = 0.2. q' = 0.2. The left, middle, and right solu- 

tions are plotted. 



3. Continuous limit and probability distribution function 

In this section, we consider the continuous limit of eq.(2.7). It is required to take the 
limit N — > oo with non zero correlation because the probability distribution function of the 
uncorrelated variables is a binomial distribution. Its continuous limit reduces to a trivial delta 
function. We need to consider the continuous limit with fixed Pd and p. Writing explicitly, 

p = (< SiSj > — < Si >< Sj >)/Pd(l — Pd) is calculated as 

p = {p 2 [(i+( 9 ' 2 -2 9 ')(i-p)r-(i-v(i- P )+g' 2 (i-p) 2 r] 

+ (1 - p)\l + (q 2 - 2q)p) N - 2 - (1 - 2qp + qV)^ 1 ] 
-2p(l - p) [(1 - p)q(l - pq) N - 2 - pq'(l - q'{l - p)) N ~ 2 
+ (l- q ')(l- q )(l-qp- q '(l-p)) N ^ 

_(1 _ g '(l -p)) N ~\l - qpf^W/Pdil - P d ), (3.1) 

where 

Pd = p(l - q'(l - p)) N ^ + (1 -p)(l - (1 - qp) N ^). (3.2) 

There are three terms in eq.(3.1). The first term is derived from < XiXj >, which is propor- 
tional to p 2 . The second term is derived from < (1 — Xi)(l — Xj) >, which is proportional 
to (1 — p) 2 . The last term is derived from < Xi(l — Xj) > and < (1 — Xj)Xj >, which is 
proportional to 2p{\ — p). At least one term must be non-zero in the continuous limit in order 
to retain the correlation. In order to fix Pd in the limit N — > oo, it is necessary to set the 
parameters p x q or (1 — p) x q' to be proportional to 1/N due to the presence of the iVth power 
in eq.(3.2). To satisfy this condition, we must set p,q,q' such that the non zero correlation 
is maintained. Using a proportional coefficient a, if we set p = a/N, the resulting expression 
corresponds to the left solution in the previous section. The correlation is maintained due to 
the first term of eq.(3.1). For p = 1 — a/N, the resulting expression corresponds to the right 
solution; the second term of eq.(3.1) remains. Instead, if we set q, q' tx 1/N and p to be finite, 
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eq.(3.1) vanishes in the limit N — > oo and the correlation disappears. The resulting expression 
corresponds to the middle solution in the previous section and the model becomes a binomial 
distribution. 

In the above limit, Pj and p can be easily estimated. We set p = a/N and substitute this 
value in eq.(3.1) and eq.(3.2). We obtain 

P d = l- e~ a \ (3.3) 



\( e <x<f _ 1) 

(3.4) 



1 - e~ a( i 



If we set p = 1 — a/N, we get 



and 



Pd = e~ aq (3.5) 



P= ~— 1 f^p 1} • (3.6) 
We see that the above two non trivial solutions can retain their non-zero correlations in the 
continuous limit. From the symmetric property of the model, we note that eq.(3.5) and eq.(3.6) 
can be derived by the substituting Pd «-> 1 — Pj, q <-> q' , and p <-> 1— p in eq.(3.3) and eq.(3.4). 

We show the profiles of the probability distributions P(k) with p = a/N and p = 1 — a/ 
N in Fig. 4 and Fig. 5. They show reflection symmetric profiles (k ^ N — /c)and have very 
singular oscillating shapes. Hereafter, we interpret the oscillating behavior based on eq.(2.7). 
The behavior of P(k) can be understood by considering each term in eq.(2.6). It is expressed 
as a summation of a^ 9 j^ (I, m). The difference between P(k) for p = a/N and p = 1 — a/N 
is only observed in the first part of eq.(2.7), p N - k ~ m +i(\ _ p^k-i+m^ p or ^ _ a /N^ 

p N- k -m + l {1 _ p)k -l +m = ^ ^ _ a j (3 ?) 

This suggests that a p ^ q ^ (l,m), which contributes significantly to P (k), should satisfy the 
condition k + m-l^N. From the second line of eq.(2.7), (1 - q 'y(k+rn-i)^ _ q ^m(N-k-m+l) ) 
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we see that I should be equal to because k — I + m ~ N. Therefore, a p ^ k with I = 
and m ~ N — k significantly contributes to P(k). The third part of eq.(2.7), (1 — (1 — 
q^N-k-m+l^k-l ^ — (1 — q'^k+m-l^N-k-m ^ ^ as non zero va } ue w jth the above condition in the 

limit N — > oo. 

We set I = and m = N — k — n (n/N -C 1) in eq.(2.7). The probability distribution 
function can be expressed as 

Jim P p=a/N (k) «jv C k x V o?=^ /N ' q ' q \l = 0,m = N-k-n), (3.8) 

where 

a P=a/N,q,q'y = Q m = N _ _ n ) = N C k X N _ k C N „ k _ n (^j (l - 

X (1 - q)( N - k - n ) n (l - (1 - q ) n ) k . (3.9) 

Instead of k, we use a normalized variable x=jfr and express a.x' q ' q (n) = ot p ^~^ N ' q ' q (I = 
0, m = N — Nx — n). The function nCnx ■ctx' q ' q (n) has a very narrow profile, and the position 
of the peak x n is given by the condition d^C^x • a"' 9,9 (n)/dx = at x = x n . We get 

77 

x n (l - q) n = (l-x n - -)(1 - (1 - g)»). (3.10) 
In the limit ./V — ► oo, we obtain the probability density function p(x) as follows. 

pW = E ^ 5 exp(-(^#), (3.11) 



where 



p(x) can be expressed as 



x n = (3.12) 
^ = ^(l-g) n (l-(l-«n. (3.13) 



P( x ) = y^PC 31 I n)P ini (n), (3.14) 

n=0 

p(x\n) = — ==expl — , (3.15 



Pim{n) = — . (3.16) 



Pini( n ) is the probability of the occurrence of n internal bad obligors and it is a Poisson 
probability function. p(x \ n) is the resulting probability density after n bad companies appear 
and infections occur from them. By the decomposition of eq.(3.14), it is easy to understand 
the oscillating behavior of p{x). In Fig. 6, p(x) is clearly decomposed into the product of the 
normal distribution p(x \ n) and Pj n j(n). 

We comment on De Finetti's theorem, which states that the joint probability function 
of TV exchangeable Bernoulli-type variables can be expressed as a mixture of the binomial 
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Fig. 6. Plot of P(k) of the left solution in Fig. 4 and the plot of the normal distribution p(x \ n)Pi n i(n) 
for < n < 20. 



distribution function ~Bi(N,p) with some mixing function /(p). 16-18 In the limit N — ► oo, 
~Bi(N,p) becomes the delta function 5(x — p) with a suitable normalization. Our model should 
be expressed by such a mixture in the N — > oo limit. In eq.(3.11), p(x \ n) becomes 5(x — x n ) 
in the limit N —> oo. The mixing function f(p) is estimated as 

f(.P) = ^2 p ini{n)5{p ~ x n )- 



We can similarly derive p(x) for the right solution [p = 1 — a/N). The result is 

(x X n ) 



p(x) 



<7 n 



E 

n=0 



exp 



2al 



p(x | n)P in i{n) 



(1 - qT, 



(3.17) 

(3.18) 

(3.19) 
(3.20) 



The probability of the occurrence of n good obligors is given by the Poisson probability 
function Pj n j(n). p(x \ n) denotes the conditional probability density function for x. 

4. Comparison with implied default distribution 

We calibrate our model based on the premiums of iTraxx-CJ (Series 2). iTraxx-CJ is an 
equally weighted portfolio of N = 50 credit default swaps (CDSs) on Japanese companies. 
The interesting point is that they are divided into several parts (called as "tranches"). The 
tranches have priorities that are defined by the attachment point ai and detachment point 
an- The protection seller agrees to cover all the losses between aLKxotai and ajjKTatal, where 
KTotai is the initial total notional of the portfolio. That is, if the loss is below aiKj-otah the 
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Fig. 7. Plot of default probability distribution for the four cases and the implied distribution on 
August 30, 2005. 



tranche does not cover it. Only when it exceeds a^Rxotah the tranche begins to cover it. If it 
exceeds anK^otah the notional becomes zero. iTraxx-CJ has five tranches and their attachment 
and detachment points are {0%,3%}, {3%, 6%},{6%, 9%}, {9%, 12%}, and {12%, 22%}. In 
addition, there is an index with {0%,100%}. We denote these tranches as {a t L ,a l H }(i = 
1, 2, • • ■ ,6) and the initial notional as Nq = (a l H — a l L ) x Kxotai- 

The premium of these tranches depends on the expectation values of the final notional 
principal < >, where < > denotes the average over the probability loss functions of 
the portfolio. If there are k defaults, the loss is expressed as the difference of the notional 
principal time and k times the recovery rate R, set at 35%, subtracted from 1. The final 
notional principal for the i— th tranche is 



N l T (k) 



Ni k < r££i 

a H N-k(l-R) rf^]<£;<r^l (4.1) 

o *>r£Ji- 

Here, \x] denotes the smallest integer greater than x. 

The premiums of the tranches Si and of the upfront Ui are determined as follows; 

rri 

+ T<N l T > Si e- rT + (N l - < N l T >)\e~ r ^ 

= {N l Q - < N l T >)e~ r ^ . (4.2) 

Here, r is the risk-free rate of interest and we set r = 1%. The left-hand side represents the 
expected payoff of the contact and the right-hand side represents the expected loss due to 
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defaults. 15 Generally, Ui = for i > 2 and s\ = 3%. 

The premiums include information about the credit market expectations regarding the 
probability loss functions. From these premiums, it is possible to infer the probability loss 
function, which is called the "implied default distribution." It describes the probability of 
k defaults of 50 Japanese companies. We denote it by Pi mp (k). There are several ways to 
infer the implied distribution; in Fig. 7, we indicate Pi mp (k) with a dotted line based on the 
maximum entropy principle. The default probability Pd is estimated as Pd = 1.65%, and the 
default correlation p is p = 6.8%. It decreases rapidly for small k; for 9 < k < 20, it is almost 
constant at ~ 0.1%. Thereafter, Pi mp (k) rapidly decays to zero. Details of the inference process 
are provided elsewhere. 11 Many probabilistic models have been proposed to date, however they 
only yield poor fits to the implied distribution. Here, we calibrate the model parameters p, q, 
and q' and study whether or not our model effectively fits Pi mp (k). 

In the calibration, we equate the default probability Pd and default correlation p of the 
model with those of the implied ones. There are three parameters p, q, and q' in the model, 
while there is only one degree of freedom. We study P(k) for the following four cases. 

(1) Default infection only (left): q' = 0.0, q = 0.054857, and p = 0.004512. 

(2) Recovery infection only (right): q = 0.0, q' = 0.421050, and p = 0.818175. 

(3) Default infection with recovery (right): q = 0.001, q' = 0.563790, and p = 0.847362. 

(4) Default infection with recovery (right): q = 0.002, q' = 0.723940, and p = 0.864563. 

The P(k) values for the above four cases are shown in Fig. 7. In the first case, the model 
exhibits only the default infection mechanism (q' = 0), and which is indicated with a solid 
line. P(k) exhibits a sharp valley structure at k = 1 and then decreases rapidly to zero at 
k ~ 20. This profile is clearly different from that of the implied one. On the other hand, 
the model with the right solution q' > q and p ~ 1, where the recovery effect dominates 
over the default infection, the bulk profiles are smooth. They are depicted by the symbols 
n\(q = 0.0), Q(q = 0.001), and A(q = 0.002). Their profiles are closer to the implied one than 
that of the infection only case. As q increases, the tail becomes short and fat. At q = 0.001 
and q = 0.002, they look similar to the implied one. We consider the infectious recovery to be 
important to describe the implied default distribution in the framework of infectious models. 

The P(k) values for the right solutions (cases 2, 3, and 4) have another peak at k = 50. The 
peak means the probability that all the 50 companies default simultaneously. The discrepancy 
from Pim P (k) is not very serious, because the inference of the default distribution from market 
quotes depends on the details of the optimization process. Instead of the entropy maximum 
principle, if we use another method, the implied distribution might have a peak at k = 50. 

The reason why the peak appears at this position in the infectious models is that we need 
to set a large value of p for obtaining the right solution. The probability that all 50 companies 
are bad is p 50 , and it is non zero; in the case, the recovery infection does not occur because 
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Premiums 


r/i So sq sa sr (Index) 


i Tr cixx- C J 


131330 008917 002850 002000 001400 002208 


IDM 1 


0.056668 0.024918 0.007940 0.001971 0.000190 0.002203 


IDM 2 


0.107641 0.013250 0.005137 0.002432 0.000809 0.002202 


IDM 3 


0.133617 0.008996 0.003695 0.002025 0.000798 0.002203 


IDM 4 


0.096937 0.014655 0.006025 0.002792 0.000815 0.002203 


Gaussian Model 


0.207827 0.004207 0.000078 0.000000 0.000000 0.002203 



Table I. Premiums for iTraxx-CJ (Series 2) on August 30, 2005, and those derived using our model. 

there are no good companies and P(50) remains. On the other hand, the probability of k < 50 
bad companies and 50 — k good companies is Pi n t(k) = soC^ -p k (l —p) 50 ~ k . In this case, 50 — k 
good companies support k bad companies, and the resulting default number is far less than 
k. Intuitively, probability Pi n t{k) for k bad companies is shifted to the left and it changes 
to probability P(k') for k! < k defaults. As q increases, in order to fix and p, we need 
to increase p and q' . The peak at k = 50 becomes higher, and the shift of P(k) to the left 
increases. As a result, P{k) moves to the left. 

In case 1 where only default infection occurs, the distribution shifts to the right in general. 
In the case of k = where there are no bad companies, the probability (1 — p) 50 is maintained. 
The default infection does not occur and P(k) has a peak at k = 0. 

We have also compared the premiums obtained using our model with the real values. The 
first row of Table 4 lists the quotes for iTraxx-CJ (Series 2) on August 30, 2005. In the second 
row, we show the premiums derived using our model with the parameters of the above four 
cases. We find that case (3) realizes the best match with the real premiums. In the table, we 
also list the premiums based on the Gaussian copula loss function, which is a standard model 
in financial engineering. 15 As the model parameter, we use the same P^ and p values as the 
implied one's. 

5. Concluding remarks and future problems 

We have generalized the infectious default model by incorporating an infectious recovery 
effect. We have explicitly obtained the default probability function P(k) for k defaults as 
a function of model parameters p, q, and q' . We have considered the continuous limit and 
obtained the probability density function p(x) for the default ratio x = ^. We have understood 
the oscillating behavior of p(x) by decomposing it, as expressed by eq.(3.14). p(x) is expressed 
as a superposition of the occurrence of n bad obligors and the following default infection. The 
former follows a Poisson distribution and the latter obeys a normal distribution. The normal 
distributions have narrow peaks of width ~ which appear in the oscillating behavior of 
p{x). We have compared the obtained P(k) with the implied one Pi mp {k) inferred from the 
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iTraxx-CJ quotes. 

By calibrating the model parameters, the profiles appear similar with regard to the bulk 
shape. However, P(k) has a peak at k = 50. We provide an intuitive explanation for it. We 
note that the principal features of our model are solvability, fitness for the implied distribution, 
and the model can be expressed as a superposition of the Poisson distributions with only three 
parameters. 

In the future, we should study the time evolution of our model. One possibility is that we 
prepare an initial configuration of Si(t = 0), whose time evolution is expressed as 

Si(t + 1) = smiuM 1 ~ Y lM l ~ S M) + (i - Si(t))(i - n j¥i (i - ^(^(t))). (5.1) 

Yij(t) and YjAt) are independent Bernoulli- type variables at each time t, and the configura- 
tion of Sj(t) is mapped to a new configuration Si(t + 1). In the original problem, the binomial 
distribution Bi(N,p) for X; is transformed into a singular oscillating P(k). We can expect 
more dynamic and complex behaviors. Furthermore, in addition to the two-body interaction 
Yij,Y-j, three-body or many body interactions might be interesting; this can be achieved by 
maintaining the integrability of the model, to determine the extent, to which such a general- 
ization is possible. In the continuous limit, the model with a continuous mixing function f(p) 
should be searched. 

The model is defined on the complete graph, where all nodes are interconnected. However, 
in recent times, the industry networks have been extensively studied and it has been shown 
that they have complex structures. 19,20 The behavior of the model on such realistic networks 
is interesting. In addition, the relation between this model and the contact process 21-23 should 
be clarified. Despite the evident similarity of our model to the contact process, the infectious 
model proposes a new approach to the description of infection. It may be that we can obtain 
the attribute of the contact process from the infectious models by considering some limit. 
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